Generalized LTB model with Inhomogeneous Isotropic Dark Energy: Observational 

Constraints 
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We consider on-center and off-center observers in an inhomogeneous, spherically symmetric, 
isocurvature (flat) concentration of dark energy with typical size of a few Gpc. Such a concen- 
tration could be produced e.g. by a recently formed global monopole with core size that approaches 
the Hubble scale. In this case we would have what may be called 'topological quintessence' in anal- 
ogy with the well-known topological inflation. We show that the minimum comoving radius romin 
of such a dark energy inhomogeneity that is consistent with the Union2 Type la supernovae (Snla) 
data at the 3a level is romin — 1.8 Gpc. As expected, the best-fit fractional dark energy density at 
the center, fix, in, approaches the corresponding ACDM value Qx,in = 0.73 for large enough values 
of the inhomogeneity radius ro (ro > 4 Gpc) . Using the Union2 data, we show that the maximum 
allowed shift r b a -max of the observer from the center of the inhomogeneity is about 0.7ro which 
respects the Copernican principle. The model naturally predicts the existence of a preferred axis 
and alignment of the low CMB multipoles. However, the constraints on r bs-max coming from the 
magnitude of the CMB dipole remain a severe challenge to the Copernican principle and lead to 
Tobs-max < 110 Mpc even for an inhomogeneity radius as large as ro = 7 Gpc. 

PACS numbers: 98.80.Es,98.65.Dx,98.62.Sb 



1. INTRODUCTION 

Detailed cosmological observations made during the 
last decade (see [l[ for some of the latest results) have 
made it clear that the minimal cosmological model 
(sCDM) based on the cosmological principle (homogene- 
ity, isotropy, validity of General Relativity) and with a 
purely radiation and matter content (baryonic and Cold 
Dark Matter (CDM)), is inconsistent with observations. 
The simplest extension of sCDM able to reconcile it with 
the observed accelerated cosmic expansion (ACDM ) in- 
volves the existence of a homogeneous, constant in time 
energy density that constitutes about 70% of the total 
energy density of the universe and is known as the cos- 
mological constant Q. A time-dependent generalization 
of this energy has been dubbed dark energy. The repul- 
sive gravitational properties of the cosmological constant 
dominate on large scales and may induce the observed 
acceleration. Due to its simplicity and its consistency 
with most cosmological observations ACDM is currently 
the standard cosmological model, replacing sCDM, which 
held that status in the early 90's, when observational 
data were far less extensive and accurate. 

However, ACDM is also faced with some challenges, 
which originate from both theoretical and observational 
considerations. The two main theoretical challenges of 
the model are fine tuning and lack of theoretical motiva- 
tion. Indeed, in order to make the model consistent with 
observations, the energy scale of the cosmological con- 
stant must be fine-tuned to an unnaturally small value, 
120 orders of magnitude smaller than the one anticipated 
from theoretical models of quantum field theory Q • 

From the observational point of view, even though 



ACDM is in agreement with the vast majority of data, 
there is a set of observations that appear to be mildly 
inconsistent with the model at the 2 — 3a level. Most of 
these observations, which we proceed to summarize, seem 
to be related to the existence of a preferred cosmological 



• Planarity and Alignment of the CMB Low 
Multipole Moments: It has been known since 
2003 H| (see [f| @ for recent reviews) that 
the quadrupole and octopole components of the 
CMB temperature perturbation maps (the largest 
scale part of those maps) are dominated by planar 
features which are unnaturally aligned with each 
other in the direction (perpendicular to the planes) 
(I, b) ~ (250°, 60°) 0rQ. The probability that such 
combination of planarity and alignment occurs in a 
random Gaussian map is less than 0.5% 5, 6]. This 
peculiar feature has been attributed to a posteri- 
ori considerations namely within a large num- 
ber of features of a map, it is statistically antici- 
pated that some of them will present large statis- 
tical fluctuations. This type of argument would 
be stronger if the large fluctuations affected ran- 
dom uncorrelated features (for example multipoles 
with I = 53 and I = 79). In that case we would 
have an a posteriori arbitrary selection of features. 
The quadrupolc-octopole moments, however, con- 
stitute all the information we have about the CMB 
temperature fluctuation maps on the largest an- 
gular scales. Thus the observed alignment could 
also be restated by saying that the large scale fea- 
tures of the CMB maps are significantly more pla- 
nar than anticipated for random Gaussian maps. 
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The fact that the preferred plane is relatively close 
to the galactic plane (which is usually masked due 
to noise, before the construction of the power spec- 
trum is made) is probably responsible (lo| for the 
relatively low values of the angular power spectrum 
components C2, C3, which is considered by some 
authors to be an independent puzzle fill ]. 

• Large Scale Velocity Bulk Flows: Recent stud- 
ies p3 - [l4| have detected a large scale flow with a 
bulk velocity in the range [400, 1000] km/s towards 
(l,b) ~ (282°, 8°) on scales 100 - 300 Mpc. This 
detection is larger than the ACDM prediction and 
inconsistent with it at the level of 2 — 3cr. 

• Large scale alignment in the QSO optical po- 
larization data: Quasar polarization vectors are 
not randomly oriented over the sky with a proba- 
bility often in excess of 99.9%. The alignment effect 
seems to be prominent along a particular axis with 
direction (l,b) = (267°, 69°) 

• Profiles of Cluster Haloes: ACDM predicts 
shallow, low concentration and density profiles, in 
contrast to some observations which point to denser 
cluster haloes [TH, [TtJ • 

As announced, three of the above four puzzles are 
large-scale effects related to preferred cosmological di- 
rections (CMB multipole alignments, QSO polarization 
alignment and large scale bulk flows); those directions, 
moreover, appear to be not far from each other, being 
approximately normal to the axis of the ecliptic poles 
(l,b) = (96°, 30°) and lying close to the ecliptic plane 
and the equinoxes. This coincidence has triggered in- 
vestigations for possible systematic effects related to the 
CMB preferred axis but no significant such effects have 
been found [f|. 

Finally, it should be mentioned that, in addition to the 
discussed large-scale effects, the ACDM model also faces 
some issues on galactic scales (missing satellites prob- 
lem p^l - l20l | and the cusp/core nature of the central den- 
sity profiles of dwarf galaxies pH|). 

In spite of those problems, the construction of gener- 
alizations of ACDM has been motivated so far mainly 
by theoretical rather than observational considerations. 
Some of the most popular ideas among those that go be- 
yond ACDM are the following: 

• Time Dependent Dark Energy |22|: Within 
this approach, dark energy retains its homogeneous 
and isotropic character, but may present a time de- 
pendence. This time dependence aims at solving 
(or at least alleviating) the previously mentioned 
fine-tuning problem, i.e. the unnaturally small 
value of the current dark energy required by the ob- 
servations. Indeed, it turns out that a dark energy 
density which is larger at early times and evolves 
via proper dynamics to its present low value is eas- 
ier to accommodate in theoretical models. A rep- 
resentative of this class of models is quintessence, 



where the role of dark energy is played by a dynam- 
ically evolving scalar field. However, for this evo- 
lution to be consistent with the observations, the 
scalar field must be assigned an extremely (fine- 
tuned) small mass, fact that makes the additional 
complexity introduced in this class of models be of 
questionable value. 

• Modified Gravity: The repulsive gravity pro- 
vided by the cosmological constant on large scales 
can also be induced by modifications of General 
Relativity in the context of a homogeneous and 
isotropic background. Several such modifications 
have been suggested, including scalar tensor grav- 
ity [23[, DGP models that involve gravitons leak- 
ing through large extra dimensions j24| . f(R) grav- 
ity [25[ etc. These models are well motivated as 
complete geometric physical theories, but they have 
more degrees of freedom than ACDM and they also 
require significant fine-tuning in order to be con- 
sistent with both the observed cosmic accelerating 
expansion and the gravitational solar system con- 
straints that prevent large local departures from 
General Relativity. 

• Inhomogeneous and Isotropic L emaitre- 
Tolman-Bondi Void Models [26U281]: Within 
this framework, an additional dark energy compo- 
nent is not needed to secure consistency with the 
cosmological data that indicate accelerating expan- 
sion. The basic idea that lies behind them is to 
consider that the increased expansion rate occurs 
locally in space rather than at recent cosmological 
times, a fact that can be achieved by assuming a 
locally-reduced matter energy density (29|. Thus, 
the observer is placed close to the center of a gi- 
ant void with dimensions of a few Gpc [30j . Even 
though this approach is free of dark energy, it is by 
no means free of fine tuning. Apart from the un- 
natural assumption of giant-size Gpc voids, which 
are very unlikely to be produced in any cosmol- 
ogy, these models require the observer to be placed 
within a very small volume at the center of the 
void (about 10~ 6 of the total volume of the void). 
A slightly off-center observer, however, will natu- 
rally experience a preferred cosmological direction 
(towards the center of the void), which may help to 
resolve some of the observational puzzles of ACDM 
discussed above. 

In the present study we consider yet an alternative ap- 
proach for the generalization of ACDM . Instead of break- 
ing the time translation invariance of dark energy, as 
in the class of quintessence models, we break its space 
translation invariance (homogeneity). We thus consider 
a spherically symmetric overdensity of dark energy with 
typical scale of a few Gpc: Inhomogeneous Dark Energy 
(IDE). This class of models adds degrees of freedom to 
ACDM in a manner that is generically different from the 
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other three classes considered above. Thus, it also has 
specific advantages as compared to them. In particular: 

• Physical Mechanism: There is a well-defined 
physical mechanism that can produce this type 
of spherically symmetric, Hubble-scale IDE. It is 
based on applying the principles of topological in- 
flation [3l| to the case of late-time acceleration. Ac- 
cording to the idea of topological inflation, the false 
vacuum energy of the core of a topological defect 
can give rise to accelerating expansion if the core 
size reaches the Hubble scale when gravity starts 
dominating the dynamics. Thus, for example, a 
recently formed global monopole with appropriate 
scale of symmetry breaking and coupling could nat- 
urally produce a Hubble-scale, spherically symmet- 
ric, isocurvature dark energy overdensity. By anal- 
ogy with topological inflation, this mechanism may 
be called topological quintessence. 

• Naturally Large Scale: Whereas matter inho- 
mogeneities (voids) of a few Gpc are very unlikely 
to be formed through large scale structure mecha- 
nisms, IDE on Hubble scales is more natural and 
easy to construct. In the IDE case (also in con- 
trast to inhomogeneous matter models) , the homo- 
geneous limit (ACDM ) is consistent with observa- 
tions and one interesting question is 'What is the 
smallest scale of inhomogeneity for which this con- 
sistency remains?'. 

• More Natural Off-center Observers (less fine 
tuning): Due to the arbitrarily large scale of IDE 
(supported by observations, as we will see), there 
is less fine tuning for the location of the observer 
with respect to the center of the spherical inho- 
mogeneity. In fact, as the inhomogeneity scale ap- 
proaches the horizon the model becomes indistin- 
guishable from ACDM and there is no observational 
constraint for the location of the observer. 

• Natural Generation of a Preferred Cosmo- 
logical Axis: As we have seen, most of the obser- 
vational challenges of ACDM hint towards the exis- 
tence of a preferred cosmological direction. Such a 
direction is naturally provided for an off-center ob- 
server in an inhomogeneous matter or dark energy 
model: the direction that connects the location of 
the observer with the center of the inhomogeneity. 
The displacement of the observer naturally leads to 
an alignment of low CMB multipole moments and 
bulk velocity flows. 

The additional degrees of freedom (with respect to 
ACDM ) introduced in this class of models are similar 
to the corresponding ones in models that break the time 
translation invariance of dark energy. We also need a 
scalar field, but now with spatial rather than temporal 
dependence. Compared to matter void models, however, 
the IDE setup requires more (and less natural) degrees 



of freedom. This is the price to pay for the reduction in 
the amount of fine tuning. 

The goal of the present study is to investigate some 
of the observational consequences and constraints of IDE 
models. We focus on the constraints coming from Snla 
data and the CMB low multipole moments. We also dis- 
cuss the physical mechanism that could give rise to IDE 
(topological quintessence) , although we postpone a com- 
plete analysis of this mechanism for a future study (32j . 

The structure of this paper is the following: In the 
next section we derive the cosmological equations, light- 
like geodesies and luminosity distance for on and off- 
center observers in the presence of IDE. We also dis- 
cuss the basics of topological quintessence as a mecha- 
nism for producing inhomogeneous dark energy and de- 
rive the angular dependence of the CMB temperature 
fluctuations induced by the shift of the observer from 
the center of the inhomogeneity. In section 3 we use 
the predicted luminosity distance to constrain the model 
parameters using the Union2 Snla data [33[. In the 
same section we compare the values of the low multi- 
pole moments (dipole, quadrupole) of CMB perturba- 
tions induced by the displacement of the observer from 
the center to the corresponding observed values and de- 
rive constraints on the parameters of the model. Finally, 
in section 4 we conclude and discuss the implications and 
possible extensions of our results. In what follows we fix 
the present age of the universe to to = 13.7 Gyr. In the 
best fit ACDM model, this would correspond to fixing 
Hq = 71km/(sMpc). We also refer to a comoving radial 
distance simply as 'distance'. 

2. COSMOLOGY WITH INHOMOGENEOUS 
DARK ENERGY 

A. Cosmological Equations 

The line element of an expanding, spherically sym- 
metric spacetime is given by the Lemaitre-Tolman-Bondi 
(LTB) metric [HQ: 

ds 2 = -dt 2 + X 2 (r,t)dr 2 + A 2 (r,t) (dO 2 + sin 2 Odip 2 ) 

(2.1) 

We consider a spherically symmetric energy-momentum 
tensor: 

T£ = diag(-p(r) - p M (r,t),p r {r),pt{r),p t {r)) (2.2) 

where pm(j, t) is the matter density and we have allowed 
for a general inhomogeneous and spherically symmetric 
static fluid with: 

(2»{J = dmg(-p(r),p r (r),p t (r),p t (r)) (2.3) 

(being pt{r) the transverse pressure). This expres- 
sion is motivated by the energy-momentum of a global 
monopole, which asymptotically is of the form j34| : 

(T mon )£ = diagH ? /r 2 , - )? /r 2 , 0, 0) (2.4) 
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with 77 the symmetry-breaking scale related to the for- 
mation of the monopole. The Einstein equations G = 



J 



SiiGT^ may be obtained from equations (|2.1[) . (|2.2[) and 
read: 



, A" 

''ax 2 



t A'X' 

''AX* 



AX A 2 



A 



\AX 



8irG(p M + p) 



(2.5) 



A' = A'— 

X 



(2.6) 



and 




A'' 



AX 2 



A AX A'X' 
A + AX + AX 3 



-8irGp r 



^ = -8irGp t . 



(2.7) 



(2. 



The primes and dots denote derivatives with respect to 
r and t, respectively. Eq. (|2.6p leads to: 



X(r,t) = C(r)A'(r,t) 



A'(r,t) 
sfl - W) 



(2.9) 



where C(r) is a constant of integration depending only 
on the coordinate r and k(r) < 1 is a function associated 
with the spatial curvature. Thus, the LTB metric (|2.1j) 
takes the form: 



ds 



-dt z 



(A'{r,t)f 



dr + A z (r, i) (dr + sin^ (9dvj z 

(2.10) 



1 - k(r) 

Taking (|2~TU| to (|2~7|) we obtain: 

i 2 + 2 A/4 + fc(r) = -8irG Pr {r)A 2 , 



(2.11) 



which implies accelerating expansion for negative pres- 
sure and proper curvature profile. It is straightforward 
to check that this equation has the following first integral: 



P_ _ F(r) 
A* ~ 



8nG . . k(r) 
Pr(r) 



A 3 3 " v ; A 2 
This last equation may also be written as 
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(2.12) 



H 2 (r,t) = H 2 (r) 
provided we define 



(2.13) 
(2.14) 



H(r,t) 



Mr,t) 
~ A(r,t) 



I 

where A (r) = A(r,t n ) (to is the present time), H (r) = 
H(r 7 t ), ilx(r) = —8irGp r (r)/2>H 2 (r) (which is positive 
for negative radial pressure) and f2 c (r) = 1 — Qx(r) — 
Q,M{r). Note that the profile of the inhomogeneous dark 
energy pressure, Qx(r), is in principle arbitrary, being 
determined by the physical mechanism that induced the 
inhomogeneous overdensity. 

The conservation of the energy- momentum tensor 1 
leads to the equations: 



A' 



Pr 



,2A 



2A' 

—£-(jPr-Pt) = 0. 



(2.16) 
(2.17) 



In the special case of isotropic pressure (p r = pt =/= 0) and 
A(r, t) = a(t)r we obtain the homogeneous dark energy 
model which has been well studied. In what follows we 
will assume the presence of anisotropic pressure. 



B. A Physical Model: Topological Quintessence 
from Global Monopoles 

A physical model that can give rise to an energy- 
momentum tensor of the form (|2.3[) is a global monopole 
with a Hubble-scale core embedded in an expanding 
spacetime. Such a system is described by the action: 



-9 



' pi 



16tt 



(2.18) 



where the Higgs field has three components $° (a 
1, 2, 3), and the potential is 



F(r)EEH 2 (r)n M (r)Al(r) 



(2.15) 



V(<f>) = -A($ 2 -?7 2 ) 2 , $ = V$ a $ a 



(2.19) 
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with rj the vacuum expectation value and A a coupling 
constant. In this model, the global monopole field con- 
figuration 

$ a = $(<, r)f a = $(t, r) (sin 9 cos tp , sin 8 sin tp, cos 9) 

(2.20) 

is topologically stable and cannot decay to the true vac- 
uum $ — r\ in a continuous manner (34[. These field 
configurations are expected to form during cosmological 
phase transitions when the vacuum manifold of a field 
theory has the topology of the sphere S 2 via the Kib- 
ble mechanism [351 ] . They are predicted to form initial 
isocurvature perturbations |36| . The monopole solution 
has localized energy in a core where $ — >• with thickness 
d ~ A -1 / 2 ?/ -1 and 4> reaches asymptotically the vacuum 
expectation value n as a power law. The asymptotic form 
of the energy- momentum tensor is given by Eq. (|2.4[) . A 
recently formed global monopole would naturally have a 
core size comparable to the present Hubble scale 



A 



-1/2 



V 



(KT^eV) 



(2.21) 



For A = 0(1) we obtain rj ~ 10~ 33 eF, but this en- 
ergy scale can increase significantly if we allow for small 
enough values of the coupling constant A. Such a field 
configuration could also be considered as a mechanism 
leading to spatial variation of fundamental constants (like 
the fine structure constant) through coupling with the 
spatially varying scalar field of the monopole. In fact 
there have been recent claims about detection [3?} of 
spatial variation of the fine structure constant along a 
particular direction in the sky. 

By varying the action (|2.18p with respect to the met- 
ric and with respect to the scalar field, we obtain the 
following dynamical equations: 



8tt 

T 

m Pl 



dV(&) 
<9$ a 



(2.22) 



(2.23) 



where 



V = d^'W - 9fiu [-(d a $ a ) 2 + V(*)\ , (2.24) 

is the energy-momentum tensor of the global monopole. 
For the metric (|2.10[) , its components are 



_/$ 2 $' 2 $ 2 
" ~ P[r) ~ ~ IT + 2X(r, tf + A^£f 

T ar = 
K = Pr (r) 



y($) 



2 



4> 



12 



<l> 2 



2 2X(r,t) 2 A(r,t) 2 
Tl = Pt (r)=T}=(-- w - y2 -V W 




(2.25) 



FIG. 1. The angle ^ between the tangent vector to the 
geodesic and the unit vector along the axis connecting the 
center with the observer. 



For a static global monopole they asymptotically reduce 
to the form (|2.4[) and negative pressure is naturally ob- 
tained in the core where the potential energy dominates. 
Such a static solution has been shown to exist in flat 
spacetimes [34| and also in expanding spacetimes when 
the core size is smaller than the Hubble scale (Fig. 2a 
of Ref. [l^). Monopolcs with cores larger than a crit- 
ical fraction (of 0(1)) of the Hubble scale have been 
shown [3l], HH to become non-static, as their core starts 
expanding along with the background Hubble expansion. 
Such an expansion may also be induced by the false vac- 
uum energy of the monopole core leading to topological 
inflation. Even though there has been no study of the 
embedded monopole in an expanding background where 
matter is also present it is conceivable that a similar be- 
havior will also be present in that case i.e. a static so- 
lution will be supported for core sizes up to about the 
Hubble scale and for a larger size, the core will start ex- 
panding along with the background especially when the 
monopole energy density of the core starts to dominate 
over the matter energy density. The behavior of such 
solutions in the presence of matter is an issue under in- 
vestigation [32I ]. 



C. General strategy 

In what follows we focus on the cosmological equation 
(|2.13[) and use it to derive observational constraints for 
observers that lie at the center of the inhomogeneity and 
also for off-center observers. To achieve this goal we fol- 
low the following steps: 
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1. We assume a dark energy profile fixM in 
Eq. (|2.13|) . characterized by two main parameters: 
the scale ro of the inhomogeneity and the magni- 
tude of fix.in at the center (fix, in = fix(f = 0)). 

2. We solve Eq. (|2.13|) assuming flatness (fijvj(r) + 
fix(0 = 1) with initial conditions taken at the 
present time to as A(r,t ) = r and evolving back- 
wards in time (39j . This initial condition is a 
gauge choice that corresponds to the usual setting 
a(to) = 1 made in homogeneous models. Our as- 
sumption for flatness is based on the fact that the 
initial perturbations produced by topological de- 
fects are isocurvature while the universe is assumed 
to be flat before the formation of the Hubble-scale 
global monopolc. 

3. Using the derived solution we solve the geodesic 
equations and find the light-like geodesies t(X), 
r(A), 9(X) for on-center and off-center observers. 
The off-center observer geodesies depend on two 
additional parameters: the distance r (, s of the ob- 
server from the center and the angle £ between the 
tangent vector to the geodesic and the unit vec- 
tor along the axis connecting the center with the 
observer |4(} ■ The value of the angle £ for each di- 
rection in the sky depends on the assumed position 
(galactic coordinates (l c , b c )) of the center of the 
inhomogeneity (i.e. £ = £(l c ,b c ), see Fig.rjJ). 

Supernovae 

4. Using the derived geodesies we obtain the lumi- 
nosity distance and other cosmological observable 
distances. These distances depend on the model 
parameters namely r Q and fix, in for the central ob- 
server and in addition to r b s and (Z c , b c ) (the galac- 



tic coordinates of the center) for the off-center ob- 
server. 

5. We fit the derived forms of the luminosity distance 
to the measured luminosity distances of the Union2 
Snla data and derive the best-fit values of the pa- 
rameters ro, fix, in, r b s and l c ,b c and their la re- 
gions. 

Cosmic Microwave Background 

6. For the derivation of the additional CMB tempera- 
ture anisotropy induced by the shift of the observer 
from the center, we use the geodesic equations to 
derive the dependence of the redshift of the last 
scattering surface on the angle £. Having found 

r obs) we find the angular dependence of the 
temperature as 



where Ti s is the observed temperature of the last 
scattering surface. This is assumed to be uniform 
since we are focusing on the additional tempera- 
ture fluctuations induced by the displacement of 
the observer. We then convert it to temperature 
fluctuations from the mean and expand in spherical 
harmonics to find the multipole coefficients. Com- 
parison with the measured values of the moments 
(particularly of the dipole) leads to constraints on 
r bs ■ 

D. Evolution of the Scale-Factor A(r,t) 

We will consider the following form for the profile of 
the inhomogeneity 



£Im(t) = fiM,out + (fijW.in — fiM,out) 
fix(r) = fix.out + (fix,in - fix,out) 



l-tanh[(r-r )/2Ar] 

f + tanh(r /2Ar) 
f - tanh[(r - r )/2Ar] 
f + tanh(r /2Ar) ' 



(2.27) 
(2.28) 



Keeping in mind the anticipated physical model of a 
global monopole, we set fiM.out = U fix.out = and 
since we are considering isocurvature inhomogeneities in 
a flat universe (fi c (r) = 0), we keep fiji/.m = 1 — fix. in- In 
our formulae we keep fi c (r) and k(r), even though these 
quantities turn out to vanish for the particular profiles 
we have considered. We also fix A?' = 0.35 Gpc, corre- 
sponding to a fairly sharp transition from 'in' to 'out'. 
We are thus left with two parameters characterizing the 
profile: the radial size ro and the density at the center 
fix,in- 



In order to solve the dynamical equation (j2.13[) we first 
determine H (r) in units of Gyr -1 , assuming a simulta- 
neous Big-Bang and, as previously mentioned, an age of 
the universe to = 13.7 Gyr. Thus we use: 

l 



to J y/flM^x- 1 +n x (r)x 2 +fi c (r) 

which can be easily derived from Eq. (|2.f 3[) . We then use 
this last equation and the initial condition A(r, to) = r 
to solve Eq. (|2.f 3|) and obtain numerically the form of 
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A(r, t) and its derivatives. 

E. Luminosity distance 

1. On-center observer 

In the case of the on-center observer, spherical symme- 
try with respect to the observer is retained and therefore 
light travels radially towards the observer. This implies 
that d6 = d<p = along the geodesic. The remaining two 
light-like geodesies for radially incoming light rays are of 
the form HH 



dt 

dz 



A'(r,t) 



(l + z)A'(r, t) 



dr c-y/l — k(r) 



dz 



(l + z)A'(r,t) 



(2.30) 



(2.31) 



where c ~ 0.3 is the velocity on light in units of Gpc/Gyr 
used for consistency in order to obtain the luminosity 
distance in Gpc. Equations (|2.30[) and (|2.31l) with initial 
conditions t(z = 0) = to and r(z = 0) = determine t(z) 
and r(z). The angular diameter and luminosity distance 
are then simply given by: 

d A {z, r , n x ,in) = A(r(z),t(z)) , (2.32) 
d L (z, r 0j n XM ) = (1 + zf d A {z) . (2.33) 

The corresponding distance modulus, which is the dif- 
ference between the apparent magnitude m and the ab- 
solute magnitude M, is given by 



/' 



M = 5 log! 



1 Gpc 



40 



(2.34) 



since d^ is obtained in Gpc. It is now straightforward 
to fit the predicted distance moduli to the observed ones 
from the Union2 dataset and derive the best-fit values of 
parameters for the on-center observer. 



2. Off-center observer 

The light-like geodesies that go through the displaced 
off-center observer (i(A), r{\), 0(X), </>(A)) (A is the affine 
parameter) , are derived by solving the geodesic equations 
in a coordinate system with origin at the center of the 
inhomogeneity [40| (see Fig. [IJ . Taking into account the 
axial symmetry around the observer, it is clear that we 
will have d<fi = along the geodesic. Moreover, due to 
the the 4- velocity identity = u^u^ = (valid 

for light-like geodesies), the time geodesic reduces to a 
first order ordinary differential equation. The geodesic 
equation for 8(A) also reduces to first order thanks to the 
spherical symmetry of the inhomogeneity (conservation 
of angular momentum). Finally, the geodesic for r(A) 



is second order, but it can be split in two first order 
equations by means of the convenient definition p = ^ . 
The resulting system consists of the following five first 
order equations [!o|, |43| : 



dt 
dX 
dr 

dX 
d6 
dX 
dz 
dX 



1 - k{r) P A 2 



cp, 
J 

C T 2 



(1+*) 



V l-fe(r 



r,2 I J 2 



1 - k(r) P + A* J 



dp 
dX 



3V =2i'p< 



r- 



l-fc(r) A 2 A' 2 

k(r)' 



l_Mr) 
+ A 3 A 1 



2 - 2k(r) 



A"_ 
~A> 



where 



A2 d6 

J = A — = const = Jo • 
dX 



(2.35) 



(2.36) 



is the conserved angular momentum (whose conservation 
emerges from the geodesic of 0(X)). The initial conditions 
are specified at the time to when the photon arrives at 
the observer's position, which is given by r = r a b s and 
8 = since without loss of generality we may assume that 
the observer is displaced along the z axis. The photon 
hits the observer at an angle £ relative to the z-axis. 
According to the metric (|2.10[) . the spatial components 
of the unit vector along this axis are 



s/l - k(r obs ) 

V = ATI Tl-l 1 ' ' ) 

A'(r obs ,to) 



(2.37) 



where the three components are in the r, 8 and <p di- 
rection, respectively. The unit tangent vector u l to the 
photon path at the present time to is given by 



d.X 

7/7 



dr d8 dej) 
dX 1 dX ' dX 



1 



(p,J/A 2 ,0) (2.38) 



where u = dt/dX and the first factor ensures normaliza- 
tion, giju}v? = 1 (notice that = gijU 1 ^ in view of 
the light-like geodesic condition u^u^ = 0). The minus 
sign in the above equation emerges because we can choose 
the affine parameter A such that A = when t = to and 
(zjj)o = u o = U (X = 0) = — 1 since the time decreases 
from to as A increases. 

The angle £ is then given by: 



cos£ = gijU l v° 



A'{r obs ,to) p{X = Q) 
^l-k{r obs ) uo 



(2.39) 



Notice that at a time t ^ t the angle £ does not represent 
the angle between the geodesic tangent and the z axis but 
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the angle between the geodesic tangent and the radial 
unit vector (which coincides with the z axis unit vector 
only at t = to when 6 — 0). In what follows we refer to £ 
as the value obtained from Eq. (j2.39[) at the present time 
t . 

We can now express the initial condition for p(A) in 
terms of the incident photon angle £ as 



y/l - k(r obs ) 

P° = mi 7^T C0S £ = cos £ 

A'(r obS7 t Q ) 



(2.40) 



Also using the fact that u l evaluated at t = to is a unit 
vector, we find using the metric (|2.10[) and Eq. (|2.38|) 

J = A(r obs ,to) sin£ = r obs sin£ ( 2 -41) 

Thus, once the incident angle £ and observer shift r obs 
I 



are specified, the initial conditions to the system (|2.35[) 
arc fixed as follows 



t(0) =to 

r(0) = r obs 

9(0) = 

z(0) =0 

p(0) = Po 



(2.42) 



Thus, the numerical solution of the system (|2.35j) with 
initial conditions (|2.42|) is straightforward and leads to 
the geodesies i(A,£, r o6s ), r(A,£, r obs ), 9(X^,r obs ) and 
z(X, £, r obs ). By inverting the expression z = z(X, £, r obs ) 
the affine parameter A may be replaced by the rcdshift 
in t, r and 8. 

The luminosity distance is now obtained from the 
geodesies for the case of an off-center observer as 0, [45| 



di(z,£,r obs ) = (1 + z) 



s A(r(z...),t(z...)) 4 S in 2 9(z...) 



sin £ 



A'(r(z...),t(z...)) 2 f drjz.. 
A(r(z...),t(z...)f(l-k(r(z...)))\ 9£ 



d9{z...) 



2n 



where ... =, £, r obs and the partial derivatives with respect 
to £ are obtained numerically by evaluating the geodesies 
at slightly different values of £. It is easy to see that 
for the on-center observer we have dr ( z >t£ob S -0) _ q an j 
de ( z 't£°b S -0) _ ^ ( smce 9 = £) and therefore we re-obtain 
the expression of the luminosity distance for the on center 
observer (|2.33[) . 

In order to compare the predicted luminosity distance 
with the corresponding observed one through Snla data, 



(2.43) 



I 

it is important to derive the dependence of the angle £ 
on the coordinates of each Snla. Let us assume an equa- 
torial coordinate system, where coordinates are given by 
the right ascension a and declination S, since the direc- 
tions of the Snla of the Union2 dataset are provided in 
this coordinate system 46]. Being the angle between 
the center-observer axis and the observer-Snla direction, 
£ depends on both the coordinates of the inhomogeneity 
center (a c , S c ) and on the coordinates of each Snla (a, 5). 
It may be shown |44| that the function £(a, S, a c , S c ) is 
given by 



£ = arccos [- sin(90° - S c ) sin(90° - 6) cos(a + (360° - a c )) + cos(90° - S c ) cos(90° - S) 



(2.44) 



Thus, equation (|2.44p translates the parameter £ to the 
coordinates of the center of the inhomogeneity (a c , S c ) 
and introduces the Snla equatorial coordinates (a,S). 
This allows the fit to the Snla Union2 data and the 
derivation of the best fit parameter values (r obs , a c , S c ). 



F. CMB anisotropies for an off-center Observer 

We assume isotropic temperature on the last scatter- 
ing surface and focus only on those additional tempera- 
ture anisotropies that emerge due to the shift of the ob- 
server from the center of the dark energy inhomogeneity. 
Thus the CMB temperature measured by the off-center 
observer is given by Eq. (|2.26[) . In order to determine 



I 

Zis(£,,r obs ) we proceed as follows: 

• We find the light-like geodesies for an on-center ob- 
server with fixed To and Slx,in and determine the 
time at decoupling ti 8 = t(A; s ) where A; s is deter- 
mined by solving the equation z(Xi s ) — 1100. 

• For a given off-center observer shifted by r obs and 
a given direction £ we find the light-like geodesies 
and determine A; s (r f, s ,£) by solving the equation 
tu = t(Xi s ,r obs ,^) where ti s is the decoupling time 
as determined for the on-center observer. 

• We determine z is (£, r obs ) = z(X ls , r obs , £) where 
z(Xi s ,r obs ,t;) is the redshift along the geodesic as 
determined in the previous subsection. 
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In order to calculate the predicted temperature fluctua- 
tions we must also obtain the predicted mean value T of 
the CMB temperature T(£, r obs ). Notice that this value 
is not exactly equal to the temperature T; s at the last 
scattering surface. The mean temperature T is evalu- 
ated using Eq. flZ26} as H3] 



f(r a 



bs ) 



1 

4tt 



dQT(£,r obs ) 



111 
2 



sinC 



1 + Zi s (€, r obs ) 
(2.45) 

The predicted angular dependence of the temperature 
fluctuations may now be easily evaluated using Eq. (|2.45[) 
as: 

AT _ T(£,r obs )-f(r obs ) 

,f, {£,! r O b s) — rfi/ \ 



z{r obs ) - zis(£,,r obs ) 



1 + zis(£,,r obs ) 



where z(r obs ) is defined as 



z{r obs ) = 2 



sin£ 



o 1 + z (t r o bs ) 



- 1 



(2.46) 



(2.47) 



The predicted values of the multipole coefficients are ob- 
tained by expanding the temperature fluctuations into 
spherical harmonics as 



AT 



(2.48) 



where the values of the multipole moments ai m are ob- 
tained as 



a>lr, 



( 2lx r AT 

J J — (£, r obs ) Y^t) sin Zdtdcf, (2.49) 



Due to axial symmetry of the temperature perturbations, 
all the ai m will vanish, except those with m = 0. The 
observed CMB dipole aw is of the order 10~ 3 . By de- 
manding that the predicted value of a 10 does not exceed 
1CP 3 we may obtain an upper bound on r obs . 

There is a Newtonian analytical approximation that 
can be made in order to obtain an expression for r obs ) 
valid for small observer displacement r obs . The cos- 
mological setup corresponding to a displaced observer 
from the center of a spherical inhomogcneity may be re- 
placed [4(| E2 by that of an observer with a specific pe- 
culiar velocity v obs . This peculiar velocity is due to the 
increased value of the Hubble parameter Hi n inside the 
inhomogeneity compared to the corresponding value that 
would exist if the inhomogcneity was not present (H out ). 
It is given by 



Vobs 

c 



~^obs 



(2.50) 



For example for the profiles of Eqs. (|2.27j) . (|2.28j) with 
r = 'A.yiGpc, Sljf.in = 0.69 we have 



q ~ 0.07 



r bs 
1 Gpc 



(2.51) 



This peculiar velocity leads to the observation of 
Dopplcr shifted photon frequencies and to an anisotropic 
temperature profile T obs (£, r bs)- The ratio of this 
anisotropic temperature over the isotropic temperature 
T c observed by the central observer is [40[ 



T 6s 

~t7 



i 



Zls 



"-7.S 



l-^cos^ 1 + zi s (£,r obs ) 



zis(£,r obs ) 

(2.52) 

where the first equality is obtained by a Doppler shift 
analysis and z; SjC ~ 1100 is the redshift of the last scat- 
tering surface as detected by the central observer. For 
small q we find 



zi s (£,r obs ) ~ zj S)C (l - gcos^) 



(2.53) 



Thus, using Eq. (|2.53|) in Eq. (|2.26p we have an approx- 
imate analytic expression for the predicted temperature 
profile which may lead to the evaluation of the predicted 
low multipole moments. We have checked that our nu- 
merical results presented in the next section are in rough 
agreement with the predictions of this Newtonian analyt- 
ical model. This model correctly predicts the qualitative 
angular dependence of z\ s on £ (oc 1 — gcos£) but we 
found it to be off by a factor ~ 2 on the amplitude q for 
the models considered. 



3. NUMERICAL ANALYSIS: OBSERVATIONAL 
CONSTRAINTS 

In this section, we implement the methods described in 
the previous section to impose observational constraints 
on the parameters of the model. In order to do so we 
will use the Union2 Snla data and also COBE/WMAP 
results for the observed CMB low multipole moments. 



A. Constraints from the Union2 Snla data 

1. On-Center Observer 

The Union2 Snla data points are given, after the cor- 
rections have been implemented, in terms of the distance 
modulus 



fJ-obs(zi) = m obs (zi) - M. 



(3.1) 



The parameters of the theoretical model are determined 
by minimizing the quantity: 



N 



(Mobs(Zi) - Mth(^i, r , Clx,m) - VoY 



j) I 



(3.2) 

where N = 557 and er^j are the errors due to flux uncer- 
tainties, intrinsic dispersion of Snla absolute magnitudes 
and peculiar velocity dispersion. These errors are as- 
sumed to be Gaussian and uncorrclated. The theoretical 
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0.20 0.25 0.30 0.35 0.40 0.45 



^M,in-l _ ^X,in 



deed consistent with the data. What is more interesting 
is the fact that the range r < 1.8 Gpc is excluded by 
the Union2 Snla data at the 3a level. The range of al- 
lowed values of fix, in is 0.53 < flx,m < 0-79 at the 3a 
level. In Fig. [2] we also show the corresponding x 2 con- 
tours for the usual LTB model (only matter + curvature, 
^M,out = !)• Notice the reduced size of the allowed re- 
gion in the parameter space and also the reduced qual- 
ity of the fit (x 2 — 556.5) as compared to the IDE and 
ACDM models (x 2 ~ 541). 

In the remainder of the section we move to investigate 
the Snla and CMB results for an off-center observer. 
Since in that case there are three additional parameters, 
( r obs,£,(c( c ,6 c )), we will consider a discrete set of values 
for (ro,fix.in), as shown in Fig. [21 with special attention 
to the values (3.37 Gpc, 0.69). These values provide a 
good fit to the Union2 data for the on-center observer, 
while being distinct from the ACDM limit. 



FIG. 2. The lcr, 2a and 3cr contours in the parameter space 
ro — Q,M,m for an on-center observer, obtained from Eq. (|3.2[) . 
The parameter values (ro, fijr.in) = (3.37 Gpc, 0.69) provide a 
good fit to the data and will be often assumed when dealing 
with the off-center case. The corresponding \ 2 contours for 
the usual LTB model (only matter + curvature, £lM,out = 1) 
are also shown. Note that the relation Q,M,in = 1 — £lx,in is 
obviously only valid for the IDE model. 

distance modulus Hth{zi,r ,ftx,in) is given by Eq. (|2.34[) 
and we have allowed for an offset fj^ which is dealt with 
either by marginalization or by minimization (the choice 
does not affect the fit of the other parameters) [13] ■ 

In Fig. [5] we show the la, 2a and 3a contours in the 
parameter space ro — fiM,in (where fijvf.in = 1 — fix. in) 
obtained from Eq. (|3.2|) after minimization with respect 
to the nuisance parameter hq. As could have been an- 
ticipated, there is no upper observational limit on the 
size of the inhomogeneity, since the homogeneous limit 
ro — > oo corresponds to the ACDM model, which is in- 



2. Off-Center Observer 

In the case of the off-center observer, the theoreti- 
cal distance moduli depend not only on (z, J"o, fijt,in)j 
but also on (r obs , £(a c , S c , a, 6)) (see Eqs. (|2~M1) . ([2~55]l . 
(|2.44p ) . In order to show the effect of the additional pa- 
rameters on the luminosity distance and the correspond- 
ing distance modulus, we show in Fig. [3] the distance 
modulus as obtained from Eqs. (|2.34[) , (|2.43p in the di- 
rection of the center (£ = ir) and in the opposite direction 
(£ = 0), for various values of the observer shift r (, s and 
fixed (ro, fix,in) = (3.37 Gpc, 0.69). We also superpose 
the Union2 data. Clearly, values of r b s > 1 Gpc in the 
direction away from the center do not provide a good 
fit. However, in the direction towards the center, the fit 
is good even for large r a b s approaching the size of the 
inhomogeneity. 

In order to find the best-fit values for the parameters 
of the model, and assuming that (ro, fix, in) ai" e fixed, we 
need to minimize the following x 2 function: 



JV 



X 2 (r obs ,a c ,<j c ,^ ) = ^2 



(Hobs(zi,cn,Si) - (ith(zi,ai,6i,r b s ,a c ,5 c ) - /xo) 5 



(3.3) 



where (i t h is obtained by combining Eqs. (|2.34[) and 
(|2.43|) . The contour plot of X2min(^c,b c ) obtained for 
(ro, ^x, in) = (3.37 Gpc, 0.69) after conversion to galactic 
coordinates and minimization with respect to r b s and 
fXo is shown in Fig. 01 The best-fit coordinate values for 
the location of the center of the dark energy spherical 
inhomogeneity are (l c , b c ) = (51.7° ± 60°, -24.9° ± 65°) 



where the errors correspond to a rough estimate of the 
la region. The corresponding la range for the observer 
shift is < r obs < 573 Mpc. In Table U and in Fig. H we 
also show the best-fit directions towards the inhomogene- 
ity center obtained for the other values (ro, fix. in) shown 
in Fig. Even though these directions are in a similar 
region in the sky, the uncertainty is clearly very large. 
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FIG. 3. The distance modulus as obtained from Eqs. (|2.43[) , (|2.34[1 in the direction of the center (£ = n) and in the opposite 
direction (£ = 0), for various values of the observer shift r b a and fixed (ro,Ox,in) = (3.37 Gpc, 0.69). The offset fio has been 
set to the best-fit value corresponding to £ = n. 



Despite of the existence of three additional param- 
eters (r ;, s , oi c , S c ), the improvement of the fit com- 
pared to on-center case is marginal (for (?"o, fijt.in) = 
(3.37Gpc,0.69) we obtain A x 2 = -2.17 at the best- 
fit point {r obs ,l Cl b c ) = (573 Mpc, 51.7°, -24.9°)). The 
improvement with respect to the ACDM model is even 
smaller: we obtain A^ 2 — —0.3 for the same values of 
(ro, fijt.in)- We conclude that there is no statistically sig- 
nificant advantage of this model over ACDM . However, 
the direction of the mildly preferred axis is not far from 
the direction of the CMB dipole (l d ,b d ) ~ (264°, 48°) = 
(84°,— 48°) where we have allowed for identification of 
opposite directions. 

In Fig. [5]wc show a contour plot of the best-fit value of 
r b s as a function of the direction in the sky (given by the 
location of the inhomogeneity center) , in galactic coordi- 
nates. The largest best-fit value is r & s = 654 Mpc corre- 
sponding to the direction (l Ctmax , b CtVaax ) = (336°, -47°). 

The maximum value of the observer shift r t> s con- 
sistent with the Union2 data may be obtained by plot- 
ting x 2 ( r obs) after minimization with respect to a c , 8 C> hq. 
This plot is shown in Fig.[5]for the set of values (ro, fJjr,in) 
shown in Fig. [5] It shows that even values of r i, s as 
large as 0.7 ro provide an acceptable fit to the Snla data 
i.e. x 2 /d.o.f. ~ 1 (even though significantly worse than 
ACDM , which corresponds to \ 2 = 541). Notice that 



ro 


fix, in 


lc 


6c 


fobs 


Xmm 


1.8 


0.54 


129° 


-26° 


0.069 


552.4 


2.45 


0.62 


87.9° 


9.0° 


0.146 


543.3 


3.0 


0.67 


62.7° 


-10.8° 


0.358 


541.1 


4.0 


0.715 


31.5° 


-41.5° 


1.05 


540.1 


5.0 


0.725 


37.1° 


-43.7° 


1.94 


539.6 



TABLE I. Best-fit direction (l c , b c ) and observer displacement 
r bs from the center of the inhomogeneity, obtained for the 
various values of (ro,fix,m) shown in Fig. [21 



shifting the observer by r a b s results only in a small im- 
provement of the fit (Ax 2 < —2) with respect to the 
on-center case. 



Predicted CMB low multipole moments: 
Constraints from the CMB dipole 



Using Eqs. (|2.49p . (|2.46|) it is straightforward to obtain 
numerically aio(r b s ), the additional CMB fluctuations 
dipole moment induced by the shift r b s of the observer. 
This quantity is shown in Fig. [7] for various values of the 
inhomogeneity scale rp. In the same plot we show the 
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FIG. 4. Contour plot of xLmCc,^) in galactic coordinates obtained after minimization with respect to r a b s and /io. The x 2 
contours shown correspond to ro = 3.37 and Ox, in = 0.69, and A\ 2 represents the improvement in the fit with respect to the 
on-center case. The best-fit directions towards the inhomogeneity center for all the (ro,Qx,in) values in Fig. [2] are also shown. 
Finally, the red point indicates the largest best-fit value for the shift of the observer (see the text and Fig. [5]). 



180" 




180" 



FIG. 5. Contour plot in galactic coordinates of the best-fit value of r j, s as a function of the location of the inhomogeneity center, 
given by the coordinates (l c ,b c )). The largest best-fit (r ba = 654 Mpc) and absolute best-fit (r & s = 573 Mpc, corresponding 
to the minimum value of \2 i n Eq. ()3.3[) ) values for the observer's shift are also shown. 



value of the measured CMB dipole (dashed line) [48 
' AT\ 3.35xlCT 3 K 




(1.230 ±0.013) x 10" 



(2.52 ±0.03) x 10" 3 

(3.4) 



Clearly, for tq < 5 Gpc the observer is confined to be in a 
sphere which is a fine tuned small spatial fraction of the 
dark energy inhomogeneity {f(r ) = ( r ° bB -™ a * ) 3 ~ 10~ 6 ) 
which implies severe fine tuning of the model as in the 
case of LTB models with matter. This is shown in Fig. [5] 
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where f(ro) is plotted for values of ro between 1 and 7 
Gpc. As shown in that figure, the fine tuning starts to 
reduce (/(ro) increases) as the size of the inhomogeneity 
increases beyond 6 Gpc. We anticipate it to disappear 
as the size of the inhomogeneity reaches ~ 14Gpc (the 
comoving distance to the last scattering surface). 

A potentially interesting observation of this class of 
models is the natural alignment of the low CMB multi- 
poles since the model includes the existence of a preferred 
axis which connects the observer with the center of the 
inhomogeneity. The predicted CMB maps correspond- 
ing to the dipolc quadrupole and octopole are shown in 
Fig. IH1 (b-d) where we have selected the direction of the 
preferred axis to coincide with the direction of the ob- 
served CMB dipole and set r b s = 30 Mpc. To try to 
avoid complications related to the numerical precision 
of the computations (complications that become mani- 
fest mainly for low values of r & s ), we have computed 
the quadrupole (octopole) in the following way: we first 
calculated its value (using Eqs. (|2^49|) . (|2~S6]> ) for 30 val- 
ues of r b s between and 1.2 Gpc and then performed 
a quadratic (cubic) fit to those data. The values for 01,20 
and 030 used in Figs. |9] and [10] have been derived from 
those fits. Notice that the predicted magnitude of the 
quadrupole is comparable with the observed value. How- 
ever, the octopole is smaller than the observed value by 
a factor ~ 10~ 4 . Even though the alignment is apparent, 



the fact that all the ai m coefficients with m 7^ vanish 
destroys all the planar features on the plane perpendic- 
ular to the preferred axis. Such features could naturally 
appear if we had considered a lattice of spherical dark 
energy inhomogeneities. 

In Fig. [TU]we show the quadrupole coefficient 020 as a 
function of the observer shift r i, s . Let us recall that the 
angular power spectrum is defined through: 

(ai m a*, m ,) = Su>S mm iCi . (3.5) 

Each C; can thus be computed by averaging the squares 
of the corresponding 21 + 1 coefficients o; TO . However, in 
our case, due to the axial symmetry (and assuming that 
the observed anisotropy were due only to our shift from 
the inhomogeneity center) we have just: 



C, = 



21 + 1 



(3.6) 



where, in addition, we should multiply by T = 2.725 x 
10 6 fiK in order to get the C;'s in units of /uK 2 , which 
is the usual convention. The yellow shaded region in 
the plot represents the value measured by WMAP, C2 
211 ■ 5/iK 2 [H, where the uncertainty only includes the 
measurement error. This corresponds to: 



5C 2 



O20 = 



(2.725 x 10V?O 5 



11.9 ± 1. 



(3.7) 
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r obs (Mpc) 

FIG. 7. The predicted CMB dipole aio(r (, s ) for fix, in = 0.69 and various values of the inhomogeneity scale ro. The dashed 
line corresponds to the observed value (whose uncertainty is too small to be shown, cf. Eq. I3.4[) 




r (Gpc) 



FIG. 8. The spatial fraction /(ro) = ( r ° bs - ma * ) 3 where 
the observer needs to be confined in order to be consistent 
with the value of the observed CMB dipole, for (ro, fix, in) = 
(3.37 Gpc, 0.69). 



We also plot the 68% (gray) and 95% (blue) confidence 
level regions accordingto the more complete statistical 
analysis carried out in [6| . Notice that the observed value 
of the quadrupole is reproduced for approximately the 



same shift of the observer r ^ s as the one required to 
produce the observed dipole. This relatively large value 
of the quadrupole we have derived is in agreement with 
the results of Ref. [42[ but is significantly larger than 
the value obtained in Ref. poj for a purely matter void 
profile. 



4. DISCUSSION AND CONCLUSIONS 

We have considered a generalization of ACDM by al- 
lowing spatial instead of temporal dependence of the dark 
energy. An advantage of this class of models is that 
they naturally predict a preferred cosmological direction 
in agreement with the indications of some cosmological 
observations [f| [TH [Tj|. We have considered isocurva- 
ture spherically symmetric dark energy inhomogeneities 
motivated from the possible recent formation of global 
monopoles with Hubble-scale cores. We have found that 
such spherical dark energy inhomogeneities with radius 
less than about 2 Gpc are excluded at the 3a confidence 
level by the Union2 Snla data. Inhomogeneities with 
larger scales are consistent with the Snla data and with 
the CMB dipole provided that the observer remains close 
to the center of the inhomogeneity, at a distance smaller 
than a maximum distance r bs-max- For the Snla Union2 
data we find r t, s - ma x ~ 0.7 ro demanding \ 2 /d.o.f. ~ 1. 



15 






-1.7 





0.9 





AT/T (quadrupole 



AT/T (octopole) 



FIG. 9. The predicted CMB maps of additional temperature fluctuations in galactic coordinates corresponding to the full 
map (a), dipole (b) quadrupole (c) and octopole (d). We have selected the direction of the preferred axis to coincide with the 
direction of the observed CMB dipole and used (ro,flx,in) = (3.37 Gpc, 0.69), r a ba = 30Mpc. Notice the rapid increase of the 
maximum allowed r b s -max as we increase the inhomogeneity size ro. It is anticipated that as we reach the value ro — 13.7 Gpc 
(the distance to the last-scattering surface) the model becomes indistinguishable from the ACDM and therefore no constraints 
are imposed on the value of r b 3 using the CMB dipole. 



This is a factor four larger than the corresponding max- 
imum observer shift in a LTB model based on a matter 
void implying that the fine tuning of the observer location 
is somewhat reduced in the case of inhomogeneous dark 
energy. This is partly due to the increase on the size of 
the observationally allowed inhomogeneities. The corre- 
sponding bound obtained from the CMB dipole however 
is not affected significantly. We find r bs-max — 110 Mpc 
even for inhomogeneities as large as 7 Gpc. This value of 
fobs -max increases rapidly as the size of the inhomogene- 
ity approaches the comoving radial distance to the last 
scattering surface ro ~ 14 Gpc. In this limit, the model 
practically reduces to ACDM . 

Our results have demonstrated that Hubble-scale in- 
homogeneities of dark energy are consistent with CMB 
and Snla cosmological observations provided that the in- 
homogeneity scale is larger than about 2 Gpc and the 
observer is kept close to the center. Potential future ex- 
tensions of this project include the following: 

• Construct a concrete physical model that can 
give rise to such large-scale dark energy inhomo- 



geneities. Topological monopolc configurations dis- 
cussed in section 2B constitute a possible physical 
basis for such models. 

• In the context of the physical model mentioned 
above, one could try to derive predictions for other 
cosmological observations. For example, the cou- 
pling of clcctromagnetism or gravity to the inho- 
mogeneous fields of the topological monopole could 
lead to a natural spatial variation of fundamental 
constants [49l [50j which would be more prominent 
along the preferred direction that connects the ob- 
server to the center of the inhomogeneity. A similar 
coupling could lead to a spatial coherent rotation 
of the optical polarization of quasars [l5[ . 

• Even though the model does predict an alignment 
of the low CMB multipoles, this alignment does 
not affect the multipoles with m ^ 0. However, 
this is no longer true if we consider a superpo- 
sition of multiple spherical inhomogeneities with 
scales and distances between their centers of a few 
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FIG. 10. a) The quadrupole moment 020 as a function of the observer shift r bs, for different inhomogeneity sizes ro and 
^A'.in = 0.69. The yellow band represents the value measured by WMAP [lj] including just the observational uncertainty (cf. 
Eq. I|3.7[l ). The 68% and 95% confidence level regions for 020 according to the more complete statistical analysis carried out 
in [(J are also shown (gray and blue shaded regions, respectively). Finally, the vertical dashed lines indicate the maximum 
value of r obs (for each ro) compatible with the measured dipole, Eq. (|3.4[) . 



Gpc. Such a lattice would keep the existence of 
a preferred axis and possible introduce concentric 
cycles in the maps [U, [HJ while inducing at the 
same time alignment between multipolcs a\ m with 
m/0. The study of the observational predictions 
of such a lattice would be an interesting extension 
of this project. 

• Hubblc-scalc dark energy inhomogeneities with dif- 
ferent symmetries would also be potentially inter- 
esting. For example, global vortices would give rise 
to cylindrically symmetric inhomogeneities while a 
domain wall would give rise to planar symmetry. 

• It has been shown that galaxies in the inner re- 
gion of a spherical inhomogeneity are seen by an 
off-center observer to have a bulk velocity flow rel- 
ative to matter outside of the inhomogeneity. The 
derivation of the predicted profile of such veloc- 
ity flows in models with inhomogeneous dark en- 
ergy would also be an interesting extension of the 
present study. 



In conclusion, the introduction of spatial variation of 
dark energy density on Hubble scales constitutes a 
natural and straightforward generalization of ACDM 
which is supported by concrete physical mechanisms. 
The detailed study of this class of models and its 
comparison with cosmological observations constitutes 
an interesting task that has been initiated by the present 
study. 

Numerical Analysis Files: The Mathematica files 
used to produce the figures are available by e-mail 
upon request and they will also be available shortly at 
http://leandros.physics.uoi.gr/idel 
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